( x ) u -d from lattice QCD at nearly physical quark masses 
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We determine the second Mellin moment of the isovector quark parton distribution function {x) u -d 
from lattice QCD with Nf — 2 sea quark flavours, employing the non-perturbatively improved 
Wilson-Sheikholeslami-Wohlert action at a pseudoscalar mass m n — 157(6) MeV. The result is 
converted non-perturbatively to the RF-MOM scheme and then perturbatively to the MS scheme 
at a scale /i = 2 GeV. As the quark mass is reduced we find the lattice prediction to approach the 
value extracted from experiments. 
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Almost all visible matter is composed of protons and 
neutrons. Analyses of the scattering of cosmic ray parti- 
cles off nuclei or of results from fixed target and colliding 
hadron beam experiments require a quantitative under- 
standing of the partonic structure of nucleons. The theo- 
retical framework for this is known [TJ [2] since the incep- 
tion of quantum chromodynamics (QCD); see also [31 E]- 

Of particular importance for the experimental pro- 
grammes at the Large Hadron Collider are unpolarized 
parton distribution functions (PDFs). These, in the 
lightcone frame, parameterize the likelihood of a parton 
to carry the Bjorken momentum fraction s at a renor- 
malization scale li. While these PDFs have been mapped 
out very well from fits to experimental data, ideally one 
would wish to evaluate them directly from the underlying 
fundamental theory, QCD. 

The present method of choice is lattice QCD, where 
in principle all approximations can be removed and sys- 
tematic uncertainties controlled by taking the limits of 
infinite volume, of vanishing lattice spacing (a — > 0) and 
of physical quark masses. However, in this Monte Carlo 
simulation approach to QCD, the statistical errors and 
the reliability of the extrapolation to the physical point 
are limited by the power of available computers and the 
efficiency of numerical algorithms. Moreover, only Mellin 
moments of the PDFs can be accessed. Thus, present-day 
lattice simulation cannot compete in terms of precision 
with determinations of isovector unpolarized PDFs from 
fits to experimental photon-nucleon scattering data that 
have been collected over decades of dedicated effort. For 
a summary of the present status of PDF parametriza- 
tions, see [5]. 

The possibility to predict averages over the momen- 
tum fraction, however, is complementary to experimental 
measurements that can only cover a limited range of x 
values. In particular, the strangeness and gluonic PDFs 
are determined rather indirectly from experiment, with 
large uncertainties [5]. Consequently, lattice QCD input 
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is already on the verge of becoming essential to constrain 
these and other less well known quantities, e.g., the pion 
nucleon a term [TVlllj. the strangeness fraction of the 
mass of the proton fx s [8Tll3|. the strange quark con- 
tribution to the spin of the proton As + As [131 H3] or 
individual (valence and sea) quark contributions to the 
proton's momentum (x) u , (x)d and (x) s |15j . Naturally, 
for lattice predictions of such quantities to be trusted, 
lattice QCD needs to demonstrate its ability to repro- 
duce known observables within similar or smaller errors. 

We focus on the second moment of the isovector PDF 
of the proton (x) u -d, i.e., on the difference between the 
average momentum fractions carried by the up and the 
down quarks. In the isospin symmetric limit, discon- 
nected quark line contributions cancel in this quantity, 
which makes it particularly accessible to lattice simula- 
tions. Nevertheless, at present (x) u -d seems to be the 
one benchmark lattice observable that is least reliably 
determined. 

We summarize the present status in Fig. [I] where we 
plot (x) u -d in the MS scheme, at a scale li = 2 GeV, as 
a function of the squared pion mass, obtained by differ- 
ent lattice collaborations. Simulations at physically light 
pion masses are particularly expensive and have only re- 
cently become possible for calculating relatively simple 
observables, e.g., the light hadron spectrum j!6) . 

The left-most points (at the physical pion mass ~ 
138 MeV) are obtained from recent PDF parametriza- 
tions [TTHT^] . evolved to the renormalization scale (x = 
2 GeV. The data point next to it (triangle) is the 
main result of this report. The QCDSF Collabora- 
tion points [20l [21] stem from simulations employing the 
same non-perturbatively improved |22| Sheikholeslami- 
Wohlert Nf — 2 quark action at lattice scales ranging 
from a -1 
points [23 

mass action at maximal twist at a 



at a 



2.60 GeV to a~ L « 3.26 GeV. The ETMC 
were obtained using the Nf = 2 twisted 

> 2.22 GeV and 
3.52 GeV while the RBC-UKQCD Collabora- 



tion [24] simulated Nf = 2 + 1 flavours using the (approx- 
imately) chiral domain- wall action at a -1 « 1.7 GeV. We 
have omitted data where the renormalization is not ex- 
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FIG. 1. Pion mass dependence of (x) u _ d at fj, = 2 GeV 
from N f = 2 (new, QCDSF [20] [21], ETMC [23]) and TV/ = 
2 + 1 (RBC-UKQCD [H]) lattice QCD simulations, together 
with expectations from PDF parametrizations (NNPDF [17] . 
ABM [fg], MSTW QI]). 



FIG. 2. Plateaus and fitted values (grey bands) of ratios of 
renormalized three-point over two-point functions at /3 = 5.29 
(a -1 « 2.76 GeV) for k = 0.13632 (m, w 290 MeV) and 
k = 0.13640 (m„ « 157 MeV). 



actly known [35] from the figure; see [57] for recent 
reviews. All the previous data shown in the figure are 
compatible with a constant ~ 0.25, far above the ex- 
pected value ~ 0.16 and no dependence on the different 
(quark and gluon) lattice actions, all accurate to 0(a), 
on the lattice spacing, on the lattice volume, on including 
a strange sea quark or not or on the pion mass can be 
resolved within statistical errors. 

Covariant baryon chiral perturbation theory suggests 
(x) u -d to decrease as the physical point is approached, 
see, e.g. [28]. However, the on-set of this behaviour is 
not at all visible for > 250 MeV, with the possi- 
ble (but not significant) exception of the RBC-UKQCD 
point with the smallest mass value. Different sources 
of systematics in lattice extractions of various nucleon 
matrix elements have been discussed in the past. In par- 
ticular, contaminations of the ground state signal, due to 
excited state contributions, have recently gained promi- 
nence in the literature [29||32| : reducing these pollutions 
results in smaller (x) u -d predictions. This effect, though 
noticeable, however seems to be too small to dominantly 
contribute to the deviations by factors of about 1.5 of 
present lattice results from values obtained from PDF 
parametrizations. 

We report on our new result that is obtained simulat- 
ing Nf = 2 non-perturbatively improved Sheikholcslami- 
Wohlert fermions on top of the Wilson gauge action at 
[3 = 5.29 and k = 0.13640, corresponding to = 
157(6) MeV, on a volume of 48 3 x 64 lattice points. Set- 
ting the scale from the (chirally extrapolated) nucleon 
mass we obtain the inverse lattice spacing [7J a -1 = 
2.76(5) (6) GeV, where the errors are statistical and from 
the chiral extrapolation, respectively. Our spatial linear 
lattice dimension L w 3.43 fm is rather small in units 
of the pion mass: m^L — 2.74(3). However, in previous 
simulations at w 290 MeV, differences between re- 



sults from three volumes of linear sizes m^L w 2.5, 3.4 
and 4.2 could not be resolved [20], within statistical er- 
rors smaller than the present one. Hence, we expect finite 
volume effects to be much smaller than our statistical un- 
certainty. 

The second Mellin moment of the isovector quark dis- 
tribution in the proton is given by 

{x) u -d = / dx x \u{x) + u(x) — d(x) — d(x)] , (1) 
Jo 

where u(x), d(x) and u(x), d(x) denote the quark and 
antiquark PDFs, respectively. It is computed by creating 
a proton at a Euclidean time to = 0; destroying it at a 
time tf > and inserting the current 

u - u- d - d , ( 2 ) 

projected to zero spatial momentum, at an intermediate 
time t < tf. For details see, e.g., the review [33]. We use 
mass-degenerate, electrically neutral u and d quarks and 
therefore all disconnected quark line loops cancel. The 
resulting three-point function is normalized by dividing 
out the two-point function of the propagating proton. 
Keeping tf fixed and sufficiently large, for tf ^> t ^> 
one obtains a plateau in t from which the lattice matrix 
element can be extracted. This is then translated into 
the MS scheme at a scale /i = 2 GeV, using the renor- 
malization factor determined in |34j , which contains non- 
perturbative renormalization into the intermediate RT- 
MOM scheme and a subsequent pcrturbative conversion 
at three-loop accuracy into the MS scheme. 

In Fig. [2] we show the resulting plateau and fit to the 
renormalized (x)^_ d where we use tf — 15a ps 1.07 fm. 
We also compare this to our previous result obtained for 
the same lattice spacing at sa 290 MeV (k — 0.13632), 
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FIG. 3. Extrapolations of the generalized form factor 
A2g d (Q 2 ) to Q = for the three smallest mass points. The 
left-most points (Q = 0) are determined directly from the 
forward matrix element. 

on a 40 3 x 64 lattice 20] . The on-set and the quality of the 
plateau depend on the overlap of the interpolating field 
used to create and to destroy the proton with the phys- 
ical ground state. Details of our quark and gauge field 
smearing can be found in [5] . Here we set the number of 
Wuppertal smearing iterations to 400. The larger pion 
mass data [M] were generated with sources and sinks of 
a smaller (Jacobi) smearing radius and exhibit more cur- 
vature as a function of t/a, indicating that the ground 
state overlap is inferior in this case. This is also obvi- 
ous from a comparison of the respective two-point func- 
tions. Improving the smearing, varying tf in addition 
to t and employing more sophisticated fit functions may 
somewhat reduce previous values that were obtained at 
the larger pion masses. Initial tests however indicate that 
the difference in quality of the interpolators used to cre- 
ate the proton can only explain a fraction of the observed 
reduction of (x) u —d as the quark mass is decreased. 

In Fig. [3] we display results for the generalized form fac- 
tor as a function of the squared momentum transfer 
Q 2 for the present mass point and for the previous two 
smallest QCDSF pion masses [5D]. (0 = 5.4 corresponds 
to the somewhat finer lattice scale a -1 ~ 3.26 GeV.) To 
leading order in chiral perturbation theory these are ex- 
pected to extrapolate linearly in Q 2 to (x) u -d = ^20~ d (0)' 
see, e.g., [2HH33J- The results agree with those obtained 
from the forward matrix elements that are also displayed 
in Fig. [I] (left-most points). 

This report is based on a (computationally expensive) 
analysis of the first 1000 configurations of our recent 
Nf = 2 simulations at a pion mass « 157 MeV and 
at an inverse lattice spacing a" 1 2.76 GeV. Altogether, 
about 3500 configurations have been generated at these 



parameter values on the QPACE supercomputers |35j . 
Hence, we expect the present error to reduce by a fac- 
tor close to two, once the target statistics are reached. 
We find (x) u ~d to drop substantially, an effect that was 
not visible at > 250 MeV. Part of this may be at- 
tributed to some excited state pollution of the previous 
data points. A re-analysis of these is also in progress. 
Our value at 157 MeV, 

(x)™ d (2 GeV) = 0.207(16), (3) 

is still by 2.5 standard deviations larger than the cen- 
tral values of PDF parametrizations IT7HT9] . Fig. [I] 
however, suggests (x) u _<j to decrease steeply as the phys- 
ical pion mass is approached. Moreover, our error, 
that incorporates the insignificant renormalization un- 
certainty [34], is purely statistical otherwise. It does not 
reflect any systematics of the missing continuum limit 
and infinite volume extrapolations. We remark that the 
comparison between lattice QCD results and fits to ex- 
perimental data for (x) u _d is less clean than for the in- 
dividual moments (x) u and (x}d- These, however, would 
require disconnected quark line diagrams to be computed 
on the lattice. Experimentally, the neutron is usually 
bound in a deuteron or another nucleus. The EMC ef- 
fect [3J5] will affect PDFs at large x [37] and change their 
moments, relative to those of a free neutron. On the 
lattice side, isospin symmetry was assumed. However, a 
strong slope as a function of the pion mass may also in- 
dicate significant corrections to this approximation. We 
plan to investigate this further. 

We conclude that lattice simulations are likely to re- 
produce known moments of unpolarized PDFs in the near 
future. This prediction is based on our observation of 
a steep decrease of the value of (x) u -d towards small 
pion masses. Clearly, high statistics simulations at pion 
masses smaller than 200 MeV and with good control over 
excited state contributions are required. 
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